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ABSTRACT 


A  review  of  the  assumptions  made  in  the  theory  of  the  classical 
integrating  sphere  is  given.  The  validity  of  these  assumptions  as 
applied  to  a  modified  integrating  sphere  is  discussed  in  terms  of  a 
computer  analysis  and  experimental  data.  The  results  of  the  analysis 
indicate  that  a  possible  error  of  15  percent  may  be  introduced  into  the 
absolute  reflectance  determinations.  This  can  occur  when  the  clas¬ 
sical  theory  is  applied  to  a  modified  integrating  sphere  using  the 
angular-hemispherical  technique.  The  parameters  considered  in  the 
analysis  are  irradiance,  sphere  radius,  sample  reflectance,  specular 
component  of  the  sample,  angle  of  incidence,  and  detector  location. 
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NOMENCLATURE 

n 

A  Area  of  the  interior  sphere  wall,  cm 

Aj  Area  of  piece  “j”,  cm22 

D  Detector 

E  Diagonal  matrix  of  reflectances 

F  Form  factor  matrix 

Fjj  Form  factor  from  i  to  j 

F ik j  Form  factor  from  i  to  k  to  i 

Gij  Form  surface  factor 

Irradiance  on  piece  j  from  piece  l,  w/cm2 
I  Identity  matrix 

Io  Energy,  w 

ISi  Incident  energy  on  St,  w/cm2 

Is  Incident  energy  on  S2 ,  w/cm2 

I q  ^  Incident  energy  at  4>  when  sample  is  irradiated  at  9,  w/cm^ 

1^0  Incident  energy  at  <t>  when  sphere  wall  is  irradiated  at  yk, 

w/ cm2 

k  Detector  constant,  mv/w/cm2 

N  Number  of  pieces  of  sphere 

Pi  Energy  flux  leaving  piece  “i”,  w 

Pij  Energy  flux  leaving  piece  “i”  and  incident  on  piece  “j”,  w 

Qjg.  fcnergy  flux  from  piece  "i"  to  piece  “i”,  w 

R  Radius  of  sphere,  in. 

r  Radius  of  sample  disk,  in. 

S  Sample  surface 

St  Standard  sample 


vx 


V  Detector  output,  mv 

a  Specular  component 

j8  Azimuthal  angle,  deg 

f  Incident  energy 

0  Angle  of  incident  beam  on  sample,  deg 

p  Reflectance 

4>  Polar  angle  of  reflection,  deg 

tp  Polar  angle  of  incident  radiation,  deg 

Aw  Solid  angle  containing  reflected  radiation 


SUBSCRIPTS 

d  Diffuse  reflecting  sample 

i  Piece  “i” 

j  Piece  “j” 

l  Piece  "i" 

r  Reflected 

Si  Pertaining  to  side  1  of  sample  S 

s2  Pertaining  to  side  2  of  sample  S 

s  Specularly  reflecting  sample 

st  Standard  sample 

w  Pertaining  to  interior  sphere  wall 
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SECTION  I 
INTRODUCTION 


There  are  three  basic  types  of  instruments  which  are  used  to 
measure  the  reflectances  of  surfaces.  These  systems  are  the  Coblentz 
hemisphere,  the  heated  cavity  reflectometer,  and  the  integrating  sphere. 
They  were  evaluated  by  Dunkle  {Ref.  1)  in  a  recent  review  paper  on  re¬ 
flectance  measurements.  The  hemisphere  with  variations  has  several 
basic  problems  such  as  critical  locations  of  detector  and  sample,  non- 
uniform  detector  response  to  both  angle  of  incidence  and  image  location, 
and  optical  aberrations.  The  heated  cavity  reflectometer  is  not  oper¬ 
able  in  the  visible  range  unless  very  high  temperatures  are  maintained 
in  the  cavity.  This  is  caused  by  large  errors  which  are  introduced  by 
small  temperature  gradients  in  the  heated  cavity.  Of  these  systems, 
the  integrating  sphere  is  capable  of  the  most  precise  measurements 
over  the  largest  wavelength  range  if  a  few  precautions  are  taken  about 
the  choice  of  sphere  wall  materials  and  the  location  of  sample,  detector, 
and  entrance  port. 

The  integrating  sphere  has  been  used  to  determine  the  reflectance 
of  various  surfaces  for  many  years,  and  extensive  literature  exists  de¬ 
scribing  its  applications.  Several  investigators  have  analyzed  the  theory 
of  the  integrating  sphere  (Refs.  2  through  5).  They  have  considered  the 
case  of  a  perfect  sphere  where  the  area  of  the  apertures  for  illumination 
and  viewing  as  well  as  the  sample  port  are  negligible  compared  to  the 
area  of  the  sphere.  Hardy  and  Pineo  (Ref.  6)  discussed  the  errors  caused 
by  the  finite  size  of  the  apertures,  but  they  considered  only  the  photo¬ 
metric  method  and  did  not  mention  the  mathematical  techniques  used.  The 
effects  of  finite  apertures  on  the  reflectance  were  also  treated  by  Preston 
(Ref.  7);  however,  his  method  cannot  be  easily  extended  to  more  general 
cases. 

Since  the  theory  of  the  integrating  sphere  is  the  theory  of  multiple 
reflections  in  a  confined  system,  this  problem  of  interreflection  can  be 
formulated  as  an  integral  equation.  Moon  (Ref.  8)  gave  a  discussion  of 
this  considering  the  apertures  to  be  negligibly  small.  Jacquez  and 
Kuppenheim  (Ref.  9)  extended  Moon's  work  and  derived  an  integral  equa¬ 
tion  which  considers  the  aperture  sizes,  geometry  of  the  sample  and 
standard,  and  the  effects  of  a  specular  component  on  the  reflectance  meas¬ 
urements.  All  of  these  discussions  assumed  perfect  sphere,  perfectly 
diffuse  reflections  from  the  wall,  and  that  the  wall  was  uniformly 
irradiated.  Also,  the  sample  and  standard  were  located  on  the  wall  of  the 
sphere.  When  the  sample  is  on  the  wall,  the  angle  of  incidence  of  the 
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light  beam  is  governed  by  the  size  of  the  sphere  and  the  geometry  of 
the  entrance  aperture.  In  most  cases  the  angle  of  incidence  is  5  to 
15  deg  from  the  sample  normal. 

The  need  for  reflectance  measurements  (Refs.  10  and  11)  as  a 
function  of  angle  of  incidence  of  the  light  beam  has  led  to  modifications 
in  the  standard  integrating  sphere.  The  most  common  modification  is 
to  place  the  sample  in  the  center  of  the  sphere  and  to  rotate  the  sample 
to  obtain  the  desired  angle  of  incidence.  Such  systems  have  been  de¬ 
scribed  by  Edwards  (Ref.  11),  Wood  (Ref.  12),  and  Birkebak  (Ref.  13). 
This  location  of  the  sample  alters  the  distribution  of  the  reflected  radia¬ 
tion  within  the  sphere  and,  therefore,  the  sphere  walls  are  not  uniformly 
irradiated  as  is  assumed  in  the  classical  theory.  The  purpose  of  this 
report  is  to  formulate  the  general  problem  of  the  integrating  sphere  with 
the  sample  in  the  center  and  to  determine  the  effect  of  the  center  sample 
on  the  wall  irradiance  and  on  the  calculation  of  sample  reflectance. 


SECTION  II 

CLASSICAL  INTEGRATING  SPHERE  THEORY 

The  reflectance  of  a  surface  is  normally  defined  as  the  ratio  of  re¬ 
flected  to  incident  radiation.  Previous  investigators  have  shown  that 
the  reflectance  is  a  function  of  surface  material,  radiation  wavelength, 
angle  of  the  incident  radiation,  and  surface  roughness.  For  the  following 
discussion  the  surface  material,  the  wavelength,  and  the  surface  rough¬ 
ness  will  be  held  constant.  It  is  assumed  that  the  interior  walls  of  the 
spheres  are  diffuse  reflectors  with  a  reflectance,  Pw,  of  96  percent,  and 
that  the  detector  is  a  perfect  absorber. 

The  classical  or  perfect  integrating  sphere  (Fig.  1),  as  discussed 
by  Jacquez  and  Kuppenheim  (Ref.  9),  is  considered  to  be  a  sphere  with 
negligible  size  apertures  whose  walls  are  uniformly  irradiated  by  multiple 
reflections.  The  sample,  S,  and  the  standard,  St,,  are  curved  such  that 
they  continue  the  sphere  wall  and,  therefore,  no  correction  need  be 
applied  for  flat  surfaces. 

When  radiation,  I0,  is  incident  on  the  diffusely  reflecting  sample, 
the  reflected  radiation  is  Pd*o>  Pd  being  the  reflectance  of  the  sample. 

The  irradiance  of  any  unit  area  on  the  sphere  wall  from  this  reflection 
would  be  pdI0/A  .  Then  the  detector  response  from  the  first  reflection 
would  be  given  by 


Vi  =  k PdI0/A 


(1) 
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From  the  second  reflection,  the  radiation  striking  the  detector  would  be 
(PwPdI0)/A2  from  each  unit  area  or  summing  over  all  unit  areas  the  de¬ 
tector  response  would  be 


V2  =  k 


PwP d  l0 
A 


(2) 


From  the  third  reflection 
V,  =  k  PwP.<*l° — 


(3) 


Summing  over  all  reflections,  the  total  detector  response  would  be 


Vd  =  f  Vi  =  k  -^2-(l  +  Pw  +  Pi  +  Pw  +  •••) 
or  in  closed  form 


(4) 


(5) 


Applying  the  same  method  to  the  standard  sample  yields  a  detector 
response  of 


Vst  =  k 


(6) 


and  the  ratio  of  Eq.  (5)  to  Eq.  (6)  gives  the  equation  of  Jacquez  and 
Kuppenheim  (Ref.  9) 


(7) 


This  equation  has  been  used  by  many  investigators  for  calculating  Pd  ■ 
However,  if  the  sample  is  not  a  perfectly  diffuse  reflector,  this  equa¬ 
tion  is  not  accurate.  To  illustrate  this,  the  same  sphere  with  a  per¬ 
fectly  specular  reflecting  sample  is  considered.  The  detector  does 
not  see  the  first  reflection  because  of  the  specular  reflection  from  the 
sample.  Therefore,  Vj  =  0,  and  the  detector  response  from  the 
second  reflection  would  be 


V  -  t  £_w_^_s_o_ 

V2  -  k  A 
From  the  third 

2  T 

p  p  I 

V3  =  k  -■  w^s  0 


(8) 

(9) 


Summing  over  all  reflections,  the  total  detector  response  for  the 
specular  sample  would  be 


Vs  =  2  V; 
i 


p  p  I 

•  w  r  s  o 


p  p  I 

r  w  r  s  o 


N-1 

Pw  P s  K 


(10) 
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or  in  closed  form 


Ul) 


and  since  pw  =  pst 


The  detector  response  for  a  diffuse  sample  (Eq.  (5))  and  a  specular 
sample  (Eq.  (11))  are  different  by  a  factor  of  Pw  •  Jacquez  and 
Kuppenheim  (Ref.  9)  arrive  at  the  same  general  equation  through  a 
more  rigorous  mathematical  treatment.  These  equations  show  that  if 
a  sample  has  both  a  diffuse  and  a  specular  component  neither  equation 
will  give  an  accurate  value  for  the  reflectance,  the  greatest  error  being 
1  -  Pw  • 


SECTION  III 

MODIFIED  SPHERE,  ANGULAR-HEMISPHERICAL  TECHNIQUE 


The  most  common  modification  to  the  classical  integrating  sphere 
has  been  the  placement  of  the  sample  in  the  center  of  the  sphere  (Fig.  2). 
For  example,  a  case  where  the  angular-hemispherical  technique  (Fig.  2) 
is  employed  will  be  considered.  Radiation  strikes  the  test  surface  at 
some  given  angle,  6,  and  is  reflected  into  the  top  hemisphere.  For  this 
system,  where  the  sample  is  located  in  the  center  of  the  sphere,  the 
assumption  of  equal  irradiance  on  the  wall  is  not  valid.  This  can  be 
shown  by  considering  the  first  reflection  of  incident  radiation  from  the 
sample. 

If  the  sample  is  a  perfectly  diffuse  reflector,  the  distribution  of  the 
radiation  reflected  into  the  top  hemisphere  will  follow  the  cosine  law. 

Then  the  area  of  wall  normal  to  the  sample  surface  will  have  the  highest 
irradiance,  and  the  irradiance  will  decrease  as  a  cosine  function  (Fig.  2). 
The  bottom  hemisphere  being  shadowed  by  the  sample  will  not  be 
irradiated  by  the  first  reflection.  For  subsequent  reflections,  all  areas 
of  the  sphere  wall  will  be  irradiated.  However,  the  shadowing  effect  of 
the  sample  must  be  considered  for  precise  calculations.  Calculation  of 
wall  irradiance  requires  a  computer  because  each  form  factor  and  form 
surface  factor  must  be  obtained  for  each  unit  area,  and  for  the  desired 
accuracy  a  large  number  of  terms  must  be  used  in  the  series  in  Eq.  (4). 
Such  a  computer  program  has  been  written,  and  the  mathematical 
method  used  in  the  program  is  discussed  in  the  next  section. 
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The  case  of  the  perfect  specular  reflecting  sample  is  less  com¬ 
plicated  than  that  of  the  diffuse  reflector.  The  radiation  is  incident  on 
the  sample  at  some  angle  6  and  is  specularly  reflected  at  the  angle  <£, 
where  6  =  <f>  •  Only  a  small  area  is  irradiated  by  the  first  reflection, 
and  this  area  is  normally  quite  small  compared  to  the  total  sphere  wall 
area.  Neglecting  this  small  area,  the  equation  derived  for  the  specular 
sample  (Eq.  (11))  can  be  used  to  calculate  the  reflectance  of  the  center 
sample.  However,  the  shadow  effect  again  has  been  neglected. 


SECTION  IV 

MATHEMATICAL  METHOD 

Consider  a  modified  sphere  (Fig.  3)  which  is  of  radius  R  and  has  a 
diffuse  reflecting  interior  surface  with  reflectance  P w  •  The  centrally 
located  sample  is  considered  to  be  an  infinitely  thin  disk  of  radius  r, 
and  the  system  is  symmetrical  about  the  z-axis.  Divide  the  sphere  into 
N  pieces  (Fig.  3),  determined  by  planes  perpendicular  to  the  z-axis, 
and  label  them  from  1  to  N  .  Let  the  test  surface  of  the  sample,  S,,  be 
piece  N  +  1  and  the  other  side,  Sa,  be  piece  N  +  2  with  their  respective 
specular  components  being  a,  and  <h  . 

For  case  1,  assume  that  a,  =  a2  =  0  If  energy  of  flux  Pi  leaves 
piece  i  diffusely  and  uniformly  and  then,  considering  all  reflections, 
there  will  be  a  resulting  flux  Pq  on  piece  j.  Since  Pij  is  directly  pro¬ 
portional  to  Pi,  a  form  surface  factor  can  be  defined  as 

Gij  =  Pq/Pi  (12) 

which  will  depend  only  on  the  nature  of  surface  and  the  geometrical 
configuration  of  the  system. 

It  is  shown  in  Appendix  I  that  by  making  the  pieces  small  enough, 
that  G  is  approximated  by  an  infinite  series  of  matrices  involving  only 
the  form  factor  matrix,  F,  and  the  diagonal  matrix  of  the  reflectances, 
E.  It  is  also  shown  that 


G  =  F  (I  -  EFT1  (13) 

Thus,  G  may  be  calculated  by  summing  a  series  or  by  inverting  a 
matrix.  The  first  method  was  used  in  this  report  because  better 
accuracy  was  obtained. 

Let  a  collimated  beam  of  flux  Qa  from  a  source  on  piece  S.  be  in¬ 
cident  on  piece  j.  Then  flux 


Pi  _  FiQii 


(14) 
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will  leave  piece  i  diffusely  (and  it  will  be  assumed  uniformly)  so  that 
Pij  can  be  calculated  from  Eq.  (12).  The  irradiance  Hgj  incident  on 
piece  j,  resulting  from  ,  is  found  by  dividing  Pij  by  the  area,  Aj, 
of  piece  j.  This  results  in 


Hij  = 


Gij 


(15) 


For  case  2,  assume  that  of  the  flux  reflected  from  Sk  (k  =  1,2) 
that  ak  is  reflected  specularly  and  that  (1  -  ak)  is  reflected  diffusely. 
Now  consider  flux  leaving  piece  i .  Depending  on  the  location,  a  frac¬ 
tion,  f ikj  ,  of  this  flux  will  be  reflected  specularly  from  Sk  onto  piece  j. 
It  is  shown  in  Ref.  14  that 

fikj  =  «k/>n  +  kFikj  (16) 

where  Fikj  is  the  form  factor  from  piece  i  to  the  part  of  the  image  of 
piece  j  reflected  in  the  plane  of  the  disk,  which  can  be  seen  through  Sk.. 

Now  define  a  new  "form  factor  "F'  by 

Fij  =  Fij  +  f  i  i.  j  +  f;2j  (17) 

Also,  since  only  (1  -  ak)  of  the  flux  leaving  Sk,  leaves  diffusely  let 

Fn  +  k,j  =  (l“  «k)Fn  +  k,j  (18) 

In  this  case  F'  is  used  instead  of  F  when  calculating  G. 


As  before,  the  irradiation  resulting  from  ,  where  i  <  n,  is  given 
by  Eq.  (15).  However,  consider  the  radiation  from  P^,  n  +  k.  An  amount, 
pn+ k  d  ~  «k)  Q^,n  +  k  >  leave  Sk  diffusely.  The  factor  of  (1  -  ak)  has 
already  been  considered  by  altering  the  form  factors  according  to  Eq.  (18). 
Also,  an  amount,  P&  Pn  +  k  ak  Pji',  will  be  leaving  piece  St  diffusely.  In 
actuality  this  flux  will  be  leaving  an  area  the  same  size  as  the  area  of  the 
source.  It  will  be  assumed  that  it  leaves  piece  St  uniformly.  Although 
this  is  only  an  approximation  to  the  actual  condition,  if  we  consider  that 
the  flux  is  leaving  diffusely  and  will  be  diffusely  reflected  many  times, 
the  assumption  leads  to  a  good  approximation.  With  this  assumption  it 
can  be  shown  that  the  irradiation  on  piece  j  resulting  from  Qi,n  +  k  is 

pj  _  8n  +  k  Qj.n+k  ^n+k,j  +  ak  ,  j  ^  (19) 


The  method  of  this  analysis  is  based  on  Eq.  (19).  It  is  seen  from 
Eqs.  (16),  through  (19)  that  case  2  reduces  to  case  1  when  ak  =  0. 
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SECTION  V 
RESULTS 


With  the  aid  of  the  computer  program,  the  wall  irradiance  was 
calculated  as  a  function  of  Ps1  >  Ps2  >  al>  a2>  ® 'P  <  and  R/r  .  These 
quantities  are  identified  in  Fig.  4.  Also,  the  reflectance  of  each  test 
sample  was  calculated,  employing  Eqs.  (7)  and  (11a),  to  determine 
the  error  caused  by  unequal  wall  irradiance.  Experimental  data  were 
obtained  on  similar  samples  in  a  7- in.  integrating  sphere  for  com¬ 
parison  purposes.  These  results  are  presented  graphically  and  are 
discussed  in  the  next  section.  These  results  are  based  on  a  sphere 
which  has  been  divided  into  16  parts,  and  a  sufficient  number  of  re¬ 
flections  were  considered  to  give  an  accuracy  of  0.  01  percent  in  the 
wall  irradiance.  The  values  of  1^  given  are  the  averages  of  the  sphere 
piece  considered,  and  4>  is  equal  to  the  angle  which  designates  the  mid¬ 
point  of  the  piece. 


SECTION  VI 
DISCUSSION 

6.1  ANGULAR-HEMISPHERICAL  TECHNIQUE 

The  system  parameters  used  in  this  discussion  are  illustrated  in 
Fig.  4.  When  the  angular-hemispherical  technique  is  used,  the  radia¬ 
tion,  q  >  strikes  the  test  surface  S, ,  at  an  angle  0  measured  from  the 
surface  normal,  and  the  irradiance  on  the  wall  is  given  in  terms  of  the 
angle  <f>.  For  the  graphs  shown,  0  =  5.6  deg  unless  otherwise  stated. 

Consider  a  diffuse  reflecting  sample  which  has  a  reflectance  of 
5  percent  (pQ ,  P<  =  0.05;  a;  =  az  =  0).  The  irradiance  of  the  sample 
is  0.  7  in.  and  is  kept  constant  for  all  calculations.  The  irradiance  on 
the  wall  is  given  in  Fig.  5  as  a  function  of  0  for  integrating  spheres  with 
radii  of  3.  5,  5,  6,  and  7  in.  The  radius  of  the  sample  is  0.  7  in.  and  is 
kept  constant  for  all  calculations  in  this  report.  The  curves  indicate 
that  the  irradiance  decreases  as  cos  <f>  from  <f>  =  0  deg  to  <f>  =  90  deg, 
but  from  <P  =  90  deg  to  =  180  deg  the  irradiance  is  essentially  constant. 
If  each  reflection  is  considered  individually,  this  is  the  expected  result. 
Since  the  first  reflection  has  a  diffuse  distribution  (Fig.  4),  a,  the  area 
normal  to  the  surface  (<f>  =  0  deg)  would  have  the  higher  irradiance  and 
the  irradiance  would  decrease  to  zero  at  <j>  =  90  deg.  The  lower  hemi¬ 
sphere  would  not  be  irradiated  by  the  first  reflection.  However,  for 
each  subsequent  reflection  the  sphere  wall  would  be  equally  irradiated 
by  each  reflection  because  of  the  diffuse  reflecting  wall.  This  is  subject 
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to  a  correction  for  the  shadowing  effect  of  the  sample.  The  change  in 
irradiance  between  the  various  spheres  is  caused  by  the  increase  in 
actual  interior  surface  area  of  the  integrating  spheres  and  the  R/r 
ratio.  The  percentage  difference  in  irradiance  between  the  upper  and 
lower  hemispheres  is  the  same  in  each  case.  That  is,  if  the  irradi¬ 
ance  was  normalized  for  the  various  sphere  radii  one  would  obtain  a 
single  curve.  Also,  since  the  sample  was  assumed  to  be  a  perfectly 
diffuse  reflector,  the  same  curves  are  obtained  when  6  is  varied  from 
0  to  90  deg.  To  compare  the  perfectly  specular  reflecting  surface 
with  the  perfectly  diffuse,  the  dashed  curves  in  Fig.  5  show  the  wall 
irradiance  on  spheres  with  3.  5-  and  6.  0-in.  radii.  All  other  param¬ 
eters  are  the  same  except  that  now  a  =  1. 


Now  consider  the  case  where  the  radiation  strikes  the  wall  first. 

The  area  which  is  illuminated  by  the  incoming  radiation  is  determined 
by  the  angle  measured  from  the  normal  to  S2  (Fig.  4).  The  change 
in  irradiance  of  the  wall  under  these  conditions  varied  less  than  1.0  per¬ 
cent  from  <f>  =  0  deg  to  </>  =  180  deg.  This  variation  was  caused  primarily 
by  the  absorption  of  radiation  by  S2 .  This  observation  is  shown 
in  Fig.  6  for  three  sphere  sizes.  The  irradiance  on  S2  is  plotted 
versus  the  angle  When  the  area  normal  to  S2  is  irradiated,  the 
irradiance  on  S2  is  the  greatest  and  as  increases  to  90  deg  the  irradi¬ 
ance  decreases  as  a  function  of  the  cos  xfj .  This  is  caused  by  the  de¬ 
crease  in  area  of  S2  (form  surface  factor)  as  seen  by  the  radiating  area 
of  the  sphere  wall.  Thus,  less  radiation  from  the  first  reflection  is 
absorbed  by  S2,  and  a  slightly  higher  irradiance  is  observed  at  ijj  =  90  deg. 
The  dotted  line  is  the  irradiance  on  S2'  and  as  shown,  ISi  =  ISi  when  the 
radiating  surface  is  90  deg  from  the  normals  of  the  two  surfaces.  The 
circular  points  on  the  curve  for  the  3.  5 -in.  radius  sphere  are  experi¬ 
mental  points  measured  by  substituting  solar  cells  for  and  S2.  The 
data  have  been  normalized  with  respect  to  the  ISj,  =25  deg. 


The  following  equation  was  used  to  calculate  the  reflectance  of 
St  ,  pSi>  from  these  data. 

ps1  =  p w  20) 


Where  Pw  is  the  reflectance  of  the  wall;  1$^  is  the  irradiance  at  <f>  when 
the  sample  is  irradiated  at  6;  I is  the  irradiance  at  <f>  when  the  in¬ 
cident  radiation  strikes  the  wall  at  Since  S2  is  a  diffuse  reflecting 
sample,  will  be  a  constant  when  6  varies  from  0  to  90  deg,  and  it 

has  been  stated  that  ^  varied  less  than  1  percent  for  all  angles  of  <A 
and  <f>.  To  analyze  the  effects  of  the  parameters  on  reflectance  meas¬ 
urements,  all  the  calculations  are  for  6  =  5.  6  deg  and  <A  =  61.9  deg. 
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Using  Eq.  (20)  or  its  equivalent,  Eq.  (7),  and  Ig^  g,  P s,  was  cal¬ 
culated  for  each  value  of  <f>  for  spheres  of  radius  3.5,  5,  and  7  in. 

In  Fig.  7,  these  results  are  compared  to  P Sj  calculated  assuming  equal 
wall  irradiance.  The  reflectance  of  a  sample  is  usually  measured 
with  the  detector  located  at  ^  =  90  deg;  this  is  indicated  on  the  figure 
by  an  arrow.  As  is  shown  by  the  curves,  if  the  reflectance  of  a  sample 
is  measured  with  the  detector  at  the  indicated  point,  the  value  of  the 
reflectance  will  be  in  error.  Several  investigators  have  made  ref.ect- 
ance  measurements  with  such  systems,  have  compared  their  results 
with  other  published  data,  and  have  reported  good  agreement  with 
spheres  of  different  sizes.  However,  the  curves  (Fig.  7)  show  that 
with  the  detector  located  at  <£  near  90  deg,  the  spheres  would  give 
essentially  the  same  result  and  it  would  be  in  error  by  approximately 
4  percent.  The  circular  points  on  the  figure  are  experimental  data 
obtained  with  a  3.  5-in.  radius  sphere.  These  data  have  been  normal¬ 
ized  with  respect  to  4>  -  80  deg. 

In  Fig.  8,  P sx  calculated  from  Eq.  (7)  is  compared  to  that  expected 
when  the  sphere  wall  irradiance  does  not  vary  (solid  lines).  The  error 
in  the  measurement  caused  by  unequal  irradiance  is  shown  in  Fig.  9  for 
values  of  Psy  from  5.0  to  100  percent  for  various  size  spheres.  The 
values  of  0,  t )j,  and  4>  are  as  indicated  in  the  figure.  The  absolute  error 
is  shown  in  Fig.  10  as  a  function  of  R  and  P .  The  absolute  error  is 
shown  again  in  Fig.  11  to  illustrate  the  change  in  error  as  a  function  of 
R/r  for  different  values  of  PSj  . 

Another  factor  which  influences  the  irradiance  is  the  reflectance 
of  the  back  side  of  the  test  surface,  Ps2  (Fig.  12).  Although  the  general 
shape  of  the  distribution  curve  is  not  changed,  the  curve  is  shifted  as 
P s2  varies  from  0  to  1.  0  as  shown  by  Fig.  12.  If  the  reflectance  meas¬ 
urements  are  obtained  with  the  detector  at  =  90  deg,  the  greatest  error 
introduced  by  S2  is  2.  0  percent. 

In  precise  calculations,  the  effect  of  the  specular  component  of  the 
test  sample  must  also  be  considered.  One  would  not  expect  the  specular 
component,  a2,  of  S2  to  have  much  influence  on  the  wall  irradiance  since 
the  sample  is  essentially  irradiated  equally  from  all  angles,  and  the 
distribution  of  the  reflected  light  would  not  vary  much  if  the  sample  were 
diffuse  or  specular.  However,  the  specular  component,  alt  of  Si  is  of 
vital  importance  when  using  the  angular-hemispherical  technique.  Since 
the  incoming  radiation  is  incident  on  the  sample,  the  distribution  of  the 
radiation  on  the  top  hemisphere  is  dictated  by  and  the  angle  of  in¬ 
cidence,  Q.  Figure  13  shows  the  effect  of  a2  on  P s,  calculated  from  the 
classical  sphere  equations.  The  curve  for  a2  =  0  was  calculated  using 
the  equation  for  a  diffuse  sample,  Eq.  (7),  and  the  curve  for  =  1.0 
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from  Eq,  (11a).  The  curves  for  the  other  values  of  ^  were  obtained 
from  the  following  equation. 


Ps,  =  “  «t) 


^0 , 4> 

V-<£ 


(pw)  +  a. 


¥><£ 


(21) 


This  equation  implies  that  the  specular  and  diffuse  components  are 
completely  separable  and  that  no  interaction  occurs.  The  validity  of 
this  assumption  is  dependent  on  the  particular  surface  being  tested. 
Measurements  of  the  distribution  of  light  reflected  from  some  metal 
samples  have  shown  that  the  specular  component  could  be  treated 
separately;  that  is,  the  specular  component  was  a  sharp  spike  super¬ 
imposed  on  a  diffuse  or  spherical  distribution.  However,  other  sam¬ 
ples  have  shown  that  the  two  components  are  inseparable  because  the 
distribution  is  oval. 


The  absolute  error  in  the  determination  of  Ps,  caused  by  using 
the  classical  sphere  equations  for  the  calculation  are  shown  in  Fig.  14 
for  several  sphere  sizes.  The  largest  error  occurs  when  a,  =  0  and 
the  smallest  when  a,  =  1.  o.  The  curve  connecting  these  two  points  was 
determined  by  Eq.  (21)  and  is  subject  to  the  assumptions  previously 
stated. 

Several  investigators  have  stated  that  if  the  angular-hemispherical 
technique  is  used  and  the  detector  does  not  see  the  first  reflection  from 
the  sample,  Eq.  (11a)  will  give  the  correct  value  of  Ps,  .  Although  it  is 
true  that  the  value  of  Ps,  calculated  from  Eq.  (11a)  when  <f>  =  90  deg  is 
in  better  agreement  with  the  true  value  of  P s,  than  the  value  obtained 
from  Eq.  (7),  the  reason  given  for  using  the  equation  is  not  valid.  This 
is  shown  in  Fig.  15.  The  solid  line  represents  the  value  of  Ps,  obtained 
from  a  perfect  sphere  as  a  function  of ■  <f>  for  %  =  0  or  1.  0.  Curve  num¬ 
ber  3  shows  Ps,  as  a  function  of  <f>  calculated  from  Eq.  (11a)  for  a  sample, 
a,  =1.0,  in  a  modified  sphere.  Curve  number  2,  calculated  from  Eq.  (7), 
is  for  a  diffuse  sample.  As  indicated,  curve  number  1,  calculated  from 
Eq.  (11a)  with  a,  =  o,  gives  a  more  accurate  value  of  Ps,  at  <j>  =  90  deg 
but  it  is  purely  fortuitous.  This  can  be  shown  by  considering  the  total 
energy  in  the  system.  If  a  sample,  Ps,  =  0.  05,  is  irradiated  by  a  source 
with  I0  watts  then  the  energy  reflected  is  0.  05  I0  watts.  This  energy  is  dis¬ 
tributed  throughout  the  sphere  and  in  the  perfect  sphere  case,  the  total 
energy  is  proportional  to  the  area  under  the  curve  for  the  perfect  sphere. 
The  area  under  curve  number  3  is  proportional  to  the  energy  for  the 
a,  =1.0  case.  This  total  area  is  approximately  4  percent  less  than  that 
for  the  perfect  sphere.  This  is  to  be  expected  since  the  detector  does 
not  see  the  first  reflection  from  the  sample,  and  its  reading  would  be  re¬ 
duced  by  Pw  •  The  area  under  the  second  curve  is  0.  4  percent  less  than 
for  the  perfect  sphere,  and  the  difference  is  attributable  to  the  absorptance 


10 


AEDC-TR- 65-271 


of  S2  .  The  area  under  curve  number  1  is  2.  8  percent  greater  than  the 
perfect  sphere  case.  This  is  not  possible  since  it  would  require  an 
addition  of  energy  to  the  system.  Even  though  the  use  of  Eq.  (11a)  to 
calculate  P sl  gives  a  more  accurate  value  under  these  conditions,  it 
has  no  theoretical  justification. 

Table  I  shows  the  absolute  error  in  the  reflectance  measurement 
when  the  angular-hemispherical  technique  is  used.  These  data  are  for 
a  3.  5-in.  radius  sphere  with  samples  of  various  reflectances  as  indi¬ 
cated  in  columns  I  and  II.  The  values  listed  in  columns  III,  IV,  and  V 
were  determined  by  the  difference  in  the  known  P s2  and  the  calculated 
values.  The  calculated  values  for  column  III  were  determined  from 
Eq.  (7)  for  <f>  =  90  deg  and  6  =  5.6  deg.  The  calculated  values  used  to 
obtain  the  errors  listed  in  columns  IV  and  V  were  determined  from 
Eq.  (11a)  for  6  =  5.  6,  <f>  =  95  and  6  =  5.6,  <f>  =  130,  respectively.  The 
values  for  Psl  -  0.  96  are  experimental  points  obtained  with  a  MgO 
sample.  The  errors  shown  in  columns  III,  IV,  and  V  indicate  that  the 
classical  theory  of  the  integrating  sphere  does  not  represent  the  data 
obtained  with  a  modified  sphere.  It  is  also  apparent  from  these  errors 
that  the  location  of  the  detector  is  very  critical. 


6.2  HEMISPHERICAL-ANGULAR  TECHNIQUE 

Another  measurement  technique  which  may  be  used  with  a  modified 
sphere  is  shown  in  Fig.  16.  In  this  technique  the  radiation  enters  the 
sphere  diffusely  at  ^  =  90  deg,  and  the  sample  is  irradiated  by  reflec¬ 
tions  from  the  top  hemisphere.  The  radiation,  I;r(^,  reflected  from  the  ' 
sample,  contained  in  solid  angle  Awr  is  measured  as  a  function  of  0. 

The  radiation  incident  on  S  ,  Iha,  is  determined  by  viewing  the  radiation 
reflected  by  the  sphere  wall  into  the  same  solid  angle,  Awr  .  The  sample 
reflectance  is  then  determined  from  Eq.  (11a).  A  detailed  discussion  of 
this  technique  is  given  in  the  literature  (Refs.  11  and  13).  The  greatest 
change  in  wall  irradiance  as  <t>  =  0  ->180  deg  was  less  than  0.  3  percent 
and  occurred  when  P§1  ~  P$2  =  0.  05.  This  variance  is  caused  primarily 
by  the  shadow  effect  of  the  sample. 

A  comparison  of  absolute,  calculated,  and  measured  reflectances 
for  a  center  sample  in  a  3.  5 -in.  radius  sphere  indicated  a  variation  of 
only  0.  1  percent;  that  is,  the  data  would  follow  the  dotted  lines  shown  in 
Fig.  8.  This  implies  that  the  assumptions  and  equations  for  the  clas¬ 
sical  integrating  sphere  theory  are  valid  for  this  technique  and  that 
accurate  results  can  be  obtained  with  their  use. 
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SECTION  VII 
CONCLUSIONS 


Calculations  and  experimental  data  have  shown  that  the  assumptions 
made  in  the  classical  theory  of  the  integrating  sphere  are  not  valid  for 
an  integrating  sphere  with  a  centrally  located  sample  when  the  angular- 
hemispherical  technique  is  used.  The  error  introduced  by  unequal  wall 
irradiance  could  be  as  great  as  15  percent  depending  on  the  reflectance 
of  the  sample  and  the  location  of  the  detector.  The  parameters  P s2  •  a> 
and  the  view  factors  may  introduce  errors  of  1  to  2  percent;  the  overall 
average  error  in  absolute  measurements  would  be  about  4  percent. 

The  theory  and  assumptions  for  an  integrating  sphere,  either  clas¬ 
sical  or  modified,  were  found  to  be  valid  when  the  hemispherical - 
angular  technique  was  used.  If  the  angular-hemispherical  technique  is 
used  with  a  modified  sphere,  the  system  should  be  calibrated  with  a 
standard  sample  for  each  set  of  experimental  conditions. 
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APPENDIX  I 


From  Eq.  (12),  if  P;  is  1  then  Gi;  =  Pn  .  Therefore,  let  a  flux  of 
unit  magnitude  leave  piece  i  diffusely  and  uniformly.  Let  P-  be  the 
flux  incident  on  piece  j  which  has  been  reflected  exactly  M  times,  then, 
since  flux  of  unit  magnitude  leaving  i  is  assumed 


(I-D 


If  it  is  assumed  that  the  flux  is  redistributed  uniformly  over  a  piece 
before  being  reflected  (this  is  a  good  assumption  if  the  pieces  are 
small)  then 


M  N  +2  , 

p>i  -«?,  F< 


Clearly  the  matrix  equation  is 

F 


PM  =  t 


pM-1  EF 


M  =  0 
M  >  0 


(1-2) 


(1-3) 


Also 


00  -  w 

G  =  2  PM 

M=0 


(1-4) 


Substituting  into  Eq.  (1-4)  and  simplifying  one  obtains  Eq.  (1-3) 

G  =  F  +(  1  PMW  d-5) 

\M=0  / 

or 

G  =  F  +  G  EF  (1-6) 

Solving  this  equation  for  G  one  obtains  the  result 

G  =  F  (I  -  EF)~l  (1-7) 

Thus,  one  can  calculate  G  by  inverting  the  matrix.  A  more  complete 
discussion  of  this  method  is  given  in  Ref.  15. 
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